clear
close all;
% 椭圆生成数据
data_length = 1000;
[x,y] = oral_sample(500,300,data_length,2000,1500);
x_0 = [x,y];
[v_x,v_y] = oral_sample(100,50,data_length,0,-100);
v_0 = [v_x,v_y];
figure; hold on;
scatter(x, y, 'filled')
quiver(x, y, v_x, v_y, .8);
title('visual (m)');
% 打开csv文件
% filename = 'data.csv';
% fid1 = fopen(filename, 'w');
% 终端时间 推力约束 推力松弛变量约束下界 推力松弛变量约束上界 物理质量上下界 推力指向约束 地下飞行约束
data = cell(data_length,8);
%% 多核并行运算
parfor i = 1:data_length
    %tic;
    data(i,:) = mult_test(x_0(i,:).',v_0(i,:).');
    %time = toc;
    %print("epoch"+num2str(i)+":"+num2str(time)+"s");
    %if(i % 100)
    %print("epoch"+num2str(time)+"s");
    %end
end
%%
% fclose(fid1);